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Abstract. The NESSiE Collaboration has been setup to undertake a conclusive experiment to clarify the 
muon-neutrino disappearance measurements at short baselines in order to put severe constraints to models 
with more than the three-standard neutrinos. To this aim the current FNAL-Booster neutrino beam for a 
Short-Baseline experiment was carefully evaluated by considering the use of magnetic spectrometers at two 
sites, near and far ones. The detector locations were studied, together with the achievable performances of 
two OPERA-like spectrometers. The study was constrained by the availability of existing hardware and a 
time-schedule compatible with the undergoing project of multi-site Liquid Argon detectors at FNAL. 
The settled physics case and the kind of proposed experiment on the Booster neutrino beam would defini¬ 
tively clarify the existing tension between the iy, disappearance and the v e appearance/disappearance at 
the eV mass scale. In the context of neutrino oscillations the measurement of disappearance is a robust 
and fast approach to either reject or discover new neutrino states at the eV mass scale. We discuss an 
experimental program able to extend by more than one order of magnitude (for neutrino disappearance) 
and by almost one order of magnitude (for antineutrino disappearance) the present range of sensitivity for 
the mixing angle between standard and sterile neutrinos. These extensions are larger than those achieved 
in any other proposal presented so far. 


1 Introduction and Physics Overview 


a Now at University of Liverpool, Department of Physics, 
Oliver Lodge Laboratory, Liverpool L69 7ZE, UK 
b Also at Centre de Recherche en Astronomie Astrophysique 
et Geophysique, Alger, Algeria 


The unfolding of the physics of the neutrino is a long and 
pivotal history spanning the last 80 years. Over this period 
the interplay of theoretical hypotheses and experimental 
facts was one of the most fruitful for the progress in parti- 
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cle physics. The achievements of the last decade and a half 
brought out a coherent picture within the Standard Model 
(SM) or some minor extensions of it, namely the mixing 
of three neutrino flavour-states with three iq, V 2 and v 3 
mass eigenstates. Few years ago a non-vanishing $ 13 , the 
last still unknown mixing angle, was measured [I]. Once 
the absolute masses of neutrinos, their Majorana/Dirac 
nature and the existence and magnitude of leptonic CP 
violation be determined, the (standard) three-neutrino 
model will be beautifully settled. Still, other questions wou¬ 
ld remain open: the reason for the characteristic nature of 
neutrinos, the relation between the leptonic and hadronic 
sectors of the SM, the origin of Dark Matter and, over¬ 
all, where and how to look for Beyond Standard Model 
(BSM) physics. Neutrinos may be an excellent source of 
BSM physics and their history supports that possibility 
at length. 

There are indeed several experimental hints for devia¬ 
tions from the “coherent” neutrino oscillation picture re¬ 
called above. Many unexpected results, not corresponding 
to a discovery on a single basis, accumulated in the last 
decade and a half, bringing attention to the hypothesis of 
the existence of sterile neutrinos |2J- A White Paper [3] 
provides a comprehensive review of these issues. In partic¬ 
ular tensions in several phenomenological models grew up 
with experimental results on neutrino/antineutrino oscil¬ 
lations at Short-Baseline (SBL) and with the more recent, 
recomputed antineutrino-fluxes from nuclear reactors. 

The main source of tension originates from the ab¬ 
sence so far of any is fl disappearance signal [Jj. Limited 
experimental data are available on searches for v tl dis¬ 
appearance at SBL: the rather old CDHS experiment [5] 
and the more recent results from MiniBooNE [B], a joint 
MiniBooNE/SciBooNE analysis [7 a and the MINOS [ 8 ] 
and SuperKanriokande [9J exclusion limits reported at the 
NEUTRINO2014 conference. The tension between v e ap¬ 
pearance and v p disappearance was actually strengthened 
by the MINOS and SuperKamiokande results, even if they 
only slightly extend the disappearance exclusion re¬ 
gion set previously mainly by CDHS and at higher mass 
scale by the CCFR experiment ms. Fig. 0 shows the ex¬ 
cluded regions in the parameter space that describe SBL 

disappearance induced by a sterile neutrino. The mix¬ 
ing angle is denoted as 6 new and the squared mass dif¬ 
ference as 4\w 1 2 elu . As evident from Fig. [ll the region 
sin 2 (2 0 new ) < 0.1 is still largely unconstrained. While this 
paper was being processed by the referees of the Jour¬ 
nal new results were made available. In particular a joint 
analysis by MINOS and DAYA-BAY [llj . and the Ice- 
Cube experiment [T2]. Their results exclude part of the 
phase space sin 2 (2 8 new ) < 0.1 even if the critical region 
Am^ ew ~ 1 eV 2 is still only marginally touched, while the 
v^- v e tension from global-fits stays around 0.04-0.07 for 
sin 2 (2 9 new ) [1]. For V^ the situation is even worse as it 
will be further discussed in Section [ 6 l 

The outlined scenario promoted several proposals for 
new, exhaustive evaluations of the neutrino phenomenol¬ 
ogy at SBL. Since the end of 2012 CERN started the 
setting up of a Neutrino Platform [13j . with new infras- 



Fig. 1. (color online) The current exclusion limits on the zy, 
disappearance searches at the eV 2 scale. Blue (green) line: 
old (recent) exclusion limits on v^ from previous CDHS [5j 
and recent MiniBooNE/SciBooNE [7] measurements. The two 
filled areas correspond to the exclusion limits on the from 
CCFR [TO] and MiniBooNE-alone [ 6 ] experiments. The red 
curve corresponds to the very recent result from MINOS [5]. 
All the exclusion limits were compute at 90% C.L.. 

tructures at the North Area that, for the time being, 
does not include a neutrino beam. Meanwhile in the US, 
FNAL welcomed proposals for experiments exploiting the 
physics potentials of their two existing neutrino beams, 
the Booster and the NuMI beams, following the recom¬ 
mendations from USA HEP-P5 report 0 . Two propos¬ 
als [T51IT6] were submitted for SBL experiments at the 
Booster beam, to complement the about to start Micro- 
B 00 NE experiment They are all based on the Liquid 
Argon (LAr) technology and aim to measure the v e appea¬ 
rance at SBL, with less possibilities to investigate the // ;J 
disappearance [IS] . In this paper a complementary case 
study based on magnetic spectrometers at two different 
sites at FNAL-Booster beam is discussed, built up on the 
following considerations: 

1 . the measurement of r'JH spectrum in both normaliza¬ 
tion and shape is mandatory for a correct interpreta¬ 
tion of the v e data, even in case of a null result for the 
latter; 

2 . a decoupled measurement of v e and z/ ;l interactions 
allows to reach in the analyses the percent-level sys- 
tematics due to the different cross-sections; 

3. very massive detectors are mandatory to collect a large 
number of events thus improving the disentangling of 
systematic effects. 

This paper is organized in the following way. After the 
introduction a short overview of the NESSiE proposal is 
given. A detailed report of the studies performed on the 
constraints of the FNAL-Booster neutrino beam is drawn 
in Section [3j In Section [4] a description of the detector 
system and the corresponding outcomes are provided. The 

1 From hereafter refers to either zz M or , unless otherwise 

stated. 
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statistical analyses and the attainable exclusion-limits on 
the Vn and V M disappearances are depicted in Sections [5] 
andlm respectively. Finally conclusions are drawn. 


2 Proposal for the FNAL-Booster beam 

Assuming the use of the FNAL-Booster neutrino beam a 
detailed study of the physics case was performed along 
the lines followed when considering neutrino beams at 
CERN-PS and CERN-SPS rHTll2(7| and the approach of 
the analysis reported in [21]. A substantial difference be¬ 
tween FNAL and CERN beams is the decrease of the aver¬ 
age neutrino energy by more than a factor 2, thus making 
the study very challenging for an high Z-density detector. 
Several detector configurations were studied, investigating 
experimental aspects not fully addressed by the LAr de¬ 
tection, such as the measurements of the lepton charge on 
event-by-event basis and the lepton energy over a wide 
range. Indeed, muons from Charged Current (CC) neu¬ 
trino interactions play an important role in disentangling 
different phenomenological scenarios provided their charge 
state is determined. Also, the study of muon appearan¬ 
ce/disappearance can benefit from the large statistics of 
CC events from the primary muon neutrino beam. 

In the FNAL Booster neutrino beam the antineutrino 
contribution is rather small and it corresponds to a syste¬ 
matic effect to be taken into account. For the antineutrino 
beam the situation is rather different since a large flux of 
neutrinos is also present. From an experimental perspec¬ 
tive the possibility of an event-by-event detection of the 
primary muon charge is an added value since it allows 
to disentangle the presence/absence of new effects which 
might genuinely affect antineutrinos and neutrinos differ¬ 
ently (CP violation is possible in models with more than 
one additional sterile neutrino). This possibility is partic¬ 
ularly intriguing while running in negative horn polarity 
due to the sizable contamination from the cross-section 
enhanced interactions of parasitic neutrinos. 

The extended NESSiE proposal is available in [2?] . 
It consists in the design, construction and installation 
of two spectrometers at two sites, Near (at 110 m, on- 
axis) and Far (at 710 m, off-axis, on surface), in line with 
the FNAL-Booster beam and compatible with the pro¬ 
posed LAr detectors. Profiting of the large mass of the two 
spectrometer-systems, their performances as stand-alone 
apparatus are exploited for the ry, disappearance study. 
Besides, complementary measurements with the foreseen 
LAr-systems can be undertaken to increase their control 
of systematic errors. 

Practical constraints were assumed in order to draft 
a proposal on a conservative, manageable basis, with sus¬ 
tainable timescale and cost-wise. Well known technologies 
were considered as well as re-using of large parts of exist¬ 
ing detectors. 

The momentum and charge state measurements of mu¬ 
ons in a wide range, from few hundreds MeV/c to several 
GeV/c, over a > 50 m 2 surface, are an extremely challeng¬ 
ing task. In the following, the key features of the proposed 


experimental layout are presented. By keeping the syste¬ 
matic error at the level of 1 — 2% for the detection of the 
i/ft interactions, it will be possible to: 

• measure the v lt disappearance in a large muon-mo¬ 
mentum, range (conservatively a p^> 500 MeV/c 
cut is chosen) in order to reject existing anomalies over 
the whole expected parameter space of sterile neutrino 
oscillations at SBL; 

• collect a very large statistical sample so as to test 
the hypothesis of muon (anti)neutrino disappearance 
for values of the mixing parameter down to still un¬ 
explored regions (sin 2 (29 new ) < 0.01); 

• measure the neutrino flux at the near detector, in the 
relevant muon momentum range, in order to keep the 
systematic errors at the lowest possible values; 

• measure the sign of the muon charge to separate iy, 
from for the control of the systematic error. 

3 Beam evaluation and constraints 

For a proposal that aims to make measurements with the 
FNAL-Booster muon -neutrino beam the convolution of 
the beam features and of the muon detection constitutes 
the major constraint. An extended study was therefore 
performed. 

3.1 The Booster Neutrino Beam (BNB) 

The neutrino beam (23) is produced from protons with a 
kinetic energy of 8 GeV extracted from the Booster and 
directed to a Beryllium cylindrical target 71 cm long and 
with a 1 cm diameter. The target is surrounded by a mag¬ 
netic focusing horn pulsed with a 170 kA current at a rate 
of 5 Hz. Secondary mesons are projected into a 50 m long 
decay-pipe where they are allowed to decay in flight before 
being stopped by an absorber and the ground material. An 
additional absorber could be placed in the decay pipe at 
about 25 m from the target. This configuration, not cur¬ 
rently in use, would modify the beam properties providing 
a more point-like source at the near site and thus extra 
experimental constraints on the systematic errors. 

Neutrinos travel about horizontally at a depth of 7 m 
underground. Proton batches from the Booster contain 
~ 4.5 x 10 12 protons, have a duration of 1.6 /zs and are sub¬ 
divided into 84 bunches. Bunches are ~ 4 ns wide and are 
separated by 19 ns. The rate of batch extraction is limited 
by the horn pulsing at 5 Hz. This timing structure pro¬ 
vides a very powerful constraint to the background from 
cosmic rays. 

3.2 The Far-to-Near Ratio (FNR) 

The uncertainty on the absolute flux at MiniBooNE, 
shown in Fig. [2j top (from [23]), stays below 20% for en¬ 
ergies below 1.5 GeV, increasing drastically at larger en¬ 
ergies and also below 200 MeV. The uncertainty is dom¬ 
inated by the knowledge of proton interactions in the Be 
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Fig. 2. (color online) Top: uncertainties on the absolute flux 
of V/i. at MiniBooNE (from [23]). Bottom: uncertainties on the 
far-to-near ratio at NuMI according to different simulations 
(from [24]). 

target, which affects the angular and momentum spec¬ 
tra of neutrino parents. The result of Fig. [2] is based on 
experimental data obtained by the HARP and E910 col¬ 
laborations [23] . 

The large uncertainty on the absolute neutrino flux 
makes the use of two or more identical detectors at differ¬ 
ent baselines mandatory when searching for small disap¬ 
pearance phenomena. The ratio of the event rates at the 
far and near detectors (FNR) as function of neutrino en¬ 
ergy is a convenient variable since at first order it benefits 
from cancellation of systematics due to the common ef¬ 
fects of proton-target, neutrino cross-sections and recon¬ 
struction efficiencies. Because of these cancellations the 
uncertainty on the FNR or, equivalently, on the spectrum 
at the far site extrapolated from the spectrum at the near 
site is at the level of few percent. As an example the FNR 
for the NuMI beam is shown in bins of neutrino energy in 
Fig. i bottom (from [M]); the uncertainty ranges in the 
interval 0.5-5.0%. 

It is worth to note that, even in the absence of oscilla¬ 
tions, the energy spectra in any two detectors are different, 
thus leading to a non-flat FNR. This is especially true if 
the distance of the near detector is comparable with the 
length of the decay pipe. It is therefore necessary to mas¬ 
ter the knowledge of the FNR for physics searches. 


Assuming a transverse area for the detectors at near 
and far sites of the same order, the solid angle subtended 
by the near detector is larger than that subtended by 
the far one. Therefore, neutrinos, and mostly those from 
mesons decaying at the end of the pipe, have a higher 
probability of being detected in the near than in the far 
detector. In the far detector, on the contrary, only neutri¬ 
nos produced in a narrow forward cone are visible. This 
effect is illustrated in Fig. [3] showing the ratio of the inte¬ 
grated neutrino flux at the two locations distributed over 
the neutrino production points (radius R vs longitudinal 
coordinate Z), for a sample crossing a near (4x4 m 2 
transverse area) and a far (8x8 m 2 ) detector with front- 
surface placed at 110 and 710 m from the target, respec¬ 
tively. Neutrinos produced at large Z can be detected even 
if they are produced at relatively large angles, enhancing 
the contribution of lower energy neutrinos. On the other 
hand neutrinos from late decays come from the fast pion 
component that is more forward -boosted. The former ef¬ 
fect is the leading one so the net effect is a softer spectrum 
at the near site. 



Z (m) 


Fig. 3. (color online) Ratio between the Z-R distributions 
of neutrino production points for neutrinos crossing a near¬ 
detector transverse surface of 4 x 4 m 2 at 110 m over neutrinos 
crossing a far-detector transverse surface of 8 x 8 m 2 at 710 m. 
The Z origin is fixed at the middle of the decay pipe, R being 
the radial distance in the pipe. The near detector has a larger 
acceptance for neutrinos produced in the most downstream 
part of the decay pipe, as expected. 

In Fig. [4] (top plots) the distributions of the neutrino 
energy, E„, vs Z for neutrinos crossing the near (top left) 
and the far site (top right) are shown. The assumed detec¬ 
tor active surface is a square of 4 x 4 m 2 and 8x8 m 2 for 
the near and far detector, respectively. As anticipated, the 
energy spectrum at the near site is softer, the additional 
contribution at low energy being particularly important 
for neutrinos from late meson-decays. The distribution of 
Z is also shown in Fig. [4] for neutrinos crossing the de¬ 
tectors at the near (bottom left) and far (bottom right) 
sites. 
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Fig. 4. (color online) Distribution of E v vs Z for neutrinos 
seen in the near (top left) and far (top right) detector. Distri¬ 
bution of the Z production points for neutrinos seen in the near 
(bottom left) and far (bottom right) detectors. The assumed 
detector active surface is a square of 4 x 4 m 2 and 8x8 m 2 for 
the near and far detector, respectively. 


ments from a single detector. However, the local interac¬ 
tion is the same when two sites, near (N) and far (F), are 
used. Any estimator of the F/N ratio in terms of some 
measurable quantity correlated to the neutrino energy is 
not affected at first order by shape distortion. Therefore, 
in case of near and far detections the effect of the inter¬ 
action models becomes sub-leading. What matters is the 
convolution of the neutrino interaction model with fluxes, 
detector acceptance as well as detector composition, which 
may be different at the two sites. The amount of this sub¬ 
leading contribution depends also on the characteristic of 
the detectors: Water Cherenkov, Liquid Argon, Scintilla¬ 
tor, Iron etc. 

NESSiE is the only proposal that could plainly profit 
of its identical configuration at near and far sites (up to 
the iron composition of the corresponding slabs in the near 
and far detectors), its capability to contain the events, and 
the control of the F/N ratio in various muon momentum 
ranges with large statistics. The contribution of the in¬ 
teraction models to the systematic error of the F/N ratio 
was checked and found very small (below 1%). An example 
of the performed several checks is reported in Fig. [5] and 
Tab.[2j for two extreme cases of the axial-mass, Ma = 0.99 
and Ma = 1-35 GeV, analyzed with a full chain simulation 
(GENIE I2SJ/FLUKA plus GEANT4 applied to configu¬ 
ration 4, see next section). 


Table 1. Systematic uncertainties on the flux prediction of the 
Booster beam. 


source 

error (%) 

p- Be 7r + production 

14.7 

2 ry nucleons interaction 

2.8 

p-delivery 

2.0 

2 ry pions interaction 

1.2 

magnetic field 

2.2 

beam-line geometry 

1.0 


From these considerations it is apparent that the pre¬ 
diction of the FNR is a delicate task requiring the full 
simulation of the neutrino beam-line and the detector ac¬ 
ceptance. Moreover, the systematic uncertainties on the 
FNR parameter play a major role requiring deep investi¬ 
gation. 

The various contributions to the systematic uncertain¬ 
ties on the neutrino flux were studied in detail by the Mini- 
BooNE collaboration in (23] (Tab. [I]). At first order they 
factorize out using a double site. However, since their mag¬ 
nitude can limit the FNR accuracy, we studied in detail 
the largest contribution, which comes from the knowledge 
of the hadro-production double differential (momentum p 1 
polar angle 0) cross-sections in 8 GeV p -Be interactions. 

Other contributions are less relevant and do not prac¬ 
tically affect the FNR estimator. As an example the sys¬ 
tematic contributions due to the multi-nucleon and the 
final state interactions have been investigated. Their mod¬ 
eling can be important when measuring cross-sections or 
for the extraction of oscillation parameters with measure- 



Fig. 5. (color online) The FNR systematic distortion due 
to different modeling of the neutrino interaction. An extreme 
value has been chosen for the axial-mass, Ma = 1-35 GeV, 
compared to the standard one, Ma = 0.99 GeV. The distortion 
is shown as function of the crossed iron planes of the detector 
(see Section [4|; it corresponds to the systematic error due to 
the convolution of the neutrino interaction, the fluxes and the 
detectors’ acceptance. 


3.3 Monte Carlo beam simulation 

In order to evaluate how the hadro-production uncer¬ 
tainty affects the knowledge of the FNR in our experiment 
a new beam-line simulation was developed. The angular 
and momentum distribution of pions exiting the Be tar¬ 
get were simulated using either FLUKA (2011.2b) [26] or 
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Table 2. Systematic uncertainties on the FNR estimator as 
function of the crossed iron planes due to the convolution of the 
neutrino interaction, the fluxes and the detectors’ acceptance, 
as shown in Fig. [5] up to a muon momentum of about 1 GeV. 
Uncertainties average at 0% with a spread < 1% due to the 
discreetness of the variable used. 


nb. pi. 

(%) 

nb. pi. 

(%) 

nb. pi. 

(%) 

0 

-0.63 

4 

0.7 

8 

-0.5 

1 

0.7 

5 

0.9 

9 

-0.6 

2 

0 

6 

-0.1 

10 

-0.1 

3 

0.1 

7 

-1.1 

11 

3.7 


GEANT4 (v4.9.4 p02, QGSP 3.4 physics list). Further¬ 
more the Sanford-Wang parametrization for 7r + deter¬ 
mined from a fit of the HARP and E910 data-set in |2!T . 
was used: 


d 2 * a 

dpdfi 


cip° 2 



— At— C50(p —cgpb cos C7 8 ) 

e p b 


(1) 


with pt being the proton- beam momentum and C; (i = 
1... 7) free parameters. The additional subdominant con¬ 
tributions arising from 7r _ and K decays have been ne¬ 
glected when considering positive polarity beam configu¬ 
rations. 

For the propagation and decays of secondary mesons 
a simulation using GEANT4 libraries was developed. A 
simplified version of the beam-line geometry was adopted. 
Despite the approximations a fair agreement with the of¬ 
ficial simulation of the MiniBooNE Collaboration [23] was 
obtained. This tool is sufficient for the purpose of the 
site optimization that is described in the following. In 
order to fully take into account finite-distance effects, 
fluxes and spectra were derived after extrapolating neutri¬ 
nos down to the detector volumes without using weighing 
techniques. A total number of 7 x 10 8 protons on target 
(p.o.t.), 2.1 x 10 s p.o.t. and 1 x 10 9 pions were simulated 
with FLUKA, GEANT4 and Sanford-Wang parametriza¬ 
tion, respectively. 

In Fig. [6] the transverse distributions of neutrinos at 
a distance of 110 m from the target is shown. The root- 
mean-square (r.m.s.) of the distribution is about 5 m. The 
projected coordinate is shown in the bottom plot of Fig. [b] 
with a Gaussian fit superimposed for comparison. The plot 
indicates that a near detector placed on ground-surface 
(Y = 7 m) would severely limit the statistics (furthermore 
the angular acceptance of the far and near detectors would 
be too different). 


3.4 Choice of experimental sites 

Once the geometry and the mass of the detectors have 
been fixed additional issues affect the choice of the lo¬ 
cation of the experimental sites, near and far ones. The 
ultimate figure of merit is the power of exclusion (or dis¬ 
covery) for effects induced by sterile neutrinos in a range 
of parameters as wide as possible in a given running timer] 

2 A similar optimization process aimed to find the best lo¬ 

cation in front of the Booster beam was extensively performed 



X (m) 



Fig. 6. (color online) The neutrino beam profile at 110 m from 
the target. In the bottom plot the projection on the horizontal 
axis, A', fitted to a Gaussian curve for comparison. 


As soon as the detectors are further away from the 
target they “see” more similar spectra since the produc¬ 
tion region better approximates a point-like source. This 
helps in reducing the systematic uncertainty. On the other 
hand the larger is the distance the smaller is the size of the 
collected event sample. Moreover the lever-arm for oscil¬ 
lation studies is reduced. The reliability of the simulation 
of the neutrino spectra at the near and far sites remains 
an essential condition. This point is further addressed in 
Section 13.4.21 

On a practical basis the increasing of the depth of the 
detector sites impacts considerably on the civil engineering 
costs. Furthermore existing or proposed experimental fa¬ 
cilities (SciBooNE/LArl-ND, T150-Icarus, MiniBooNE, 
MicroBooNE, LArl, Icarus) already partially occupy pos¬ 
sible sites along the beam line p~8l . 


by the SciBooNE collaboration [27]. In that case the aim was 
either to maximize the neutrino flux or to shape out the energy 
interval for cross-section measurements. 
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Z (m) 



Fig. 7. (color online) Top: rate (a.u.) vs distance from 

the target, Z, and the depth from the ground surface. Bot¬ 
tom: mean energy (in GeV) of i/^ c interactions vs Z (m) and 
depth. The origin of the reference system is the proton-target 
upstream position. The hollow box close to the origin indicates 
the longitudinal extension of the decay pipe. 


3.4.1 Dependence of v^ c rates and energy spectra on the 
detector position 

The interaction rates and their mean energy depend 
on the distance from the proton target, as shown in Fig. [7j 
top and bottom, respectively. The horizontal axis corre¬ 
sponds to the distance (Z) from the target, the vertical 
axis to the depth from the ground surface. At a distance 
of about 700 m the rates and the mean energies are barely 
affected when moving from on-axis to off-axis positions. 
That consideration supports the possibility of placing the 
far detector on surface, thus reducing the experiment cost. 


3.4.2 Systematics in the near-to-far ratio for a set of 
detector configurations 

Six configurations were selected considering different dis¬ 
tances (110, 460 and 710 m), either on-axis or off-axis, 
and different fiducial sizes of the detectors. The config¬ 
urations’ parameters are given in Tab. [3] and illustrated 
schematically in Fig. [8] 

• Configuration 1 corresponds to two detectors on-axis 
at 110 (near) and 710 m (far) with squared active areas 


Table 3. Near- Far detectors’ configurations. Lft(f) is the dis¬ 
tance of the near (far) detector from the target. Un(f) is the 
vertical coordinate of the center of the near (far) detector with 
respect to the beam axis, which lies at about —7 m from the 
ground surface. Sjv(f) is the transverse size of the near (far) 
detector. 


config. 

L N (m) 

L f (m) 

Vn (m) 

Vf (m) 

SN (m) 

SF (ill) 

1 

110 

710 

0 

0 

4 

8 

2 

110 

710 

0 

0 

1.25 

8 

3 

110 

710 

1.4 

11 

4 

8 

4 

110 

710 

1.4 

11 

1.25 

8 

5 

460 

710 

7 

11 

4 

8 

6 

460 

710 

6.5 

10 

4 

6 
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Fig. 8. (color online) Configurations of far and near detec¬ 
tors in the Y—Z plane (see also Tab. [3]). The blue horizontal 
line marks the ground level, the vertical black lines mark the 
detectors and the red lines show the angle subtended by the 
detectors at the beginning and the end of the decay pipe. 


of 4 x 4 m 2 (near) and 8x8 m 2 (far). By selecting the 
subsample of neutrinos crossing both the near and far 
detectors the region defined in the transverse plane has 
roughly a squared shape with a significant “blurring” 
since the neutrino source is not point-like. 

• Configuration 2 make use of a reduced near detector 
area, limited to 1.25 x 1.25 m 2 , in order to increase 
the fraction of neutrinos seen both at the near and far 
sites. 

• Configurations 3 and 4 replicate the same patterns as 
1 and 2, respectively, with the far detector on surface 
and the size of the near detector defined by the off-axis 
angle (instead of being both on-axis). 

• Configurations 5 and 6 are similar to 3 and 4, respec¬ 
tively, but with the near site at a larger distance (460 
m). 

Using FLUKA, GEANT4 or the Sanford-Wang para- 
metrization for the simulation of p -Be interactions, the 
FNR was computed for each configuration (Fig. [9]). 
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In configuration 1 (on -axis detectors and a large near¬ 
detector size) the FNR increases with energy, as expected 
from the discussion in Section [3.2[ largely departing from 
a flat curve. By reducing the transverse size of the near de¬ 
tector (configuration 2) the FNR flattens out. This same 
behavior is confirmed using off-axis detectors (configura¬ 
tions 3 and 4). Even flatter FNRs are obtained in config¬ 
urations with a near detector at larger baselines (5 and 


6 ). The different behaviors are more evident in Fig. 10 


where the FNRs based on the Sanford-Wang parametri- 
zation and normalized to each other are compared. 




1 1.2 1.4 1.1 

neutrino energy (GeV) 



neutrino energy (GeV) 


neutrino energy (GeV) 


I Config. 6 I 

TTTTTTTT 



1 1.2 1.4 1.6 

neutrino energy (GeV) 


1 1.2 1.4 1.6 

neutrino energy (GeV) 


Fig. 9. (color online) Far-to-Near ratios for each consid¬ 
ered configuration. Comparison of results from the different 
hadro-productions simulations by FLUKA, GEANT4 and the 
Sanford-Wang parametrization are shown. The error bars in¬ 
dicate only the uncertainty introduced by the limited Monte 
Carlo samples of FLUKA and GEANT4. The barely visible er¬ 
ror bars on the Sanford-Wang points are due to the very large 
number of simulated pions. 


In order to estimate the impact on the FNR of the 
hadro-production uncertainties two studies were made. 

First, the difference in the hadronic models implemen¬ 
ted in the FLUKA and GEANT4 generators were looked 
at. For each configuration the FNR predictions from these 
two Monte Carlo simulations are drawn in Fig. 11 as their 
ratio. The (yellow) bands correspond to a fixed 3% er¬ 
ror on the FNR ratios between FLUKA and GEANT4. 
The two simulations agree at 1 to 3% level when an over¬ 
lapping region between the far and near detectors occurs 
(configurations 2, 4 and 6). 

Another approach was adopted to investigate the FNR 
systematic-error due to hadro-production based on the 
existing measurements and the corresponding covariance- 
error matrix of Booster Be-target replica [22]. The co¬ 
efficients Ci of the Sanford-Wang parametrization of pion 



Fig. 10. (color online) Far-to-Near ratios for the six consid¬ 
ered configurations using the Sanford-Wang parametrization. 
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Fig. 11. Ratio of the two FNRs as predicted by the FLUKA 
and GEANT4 simulations for hadro-production. The bands 
correspond to a 3% relative error on the different predictions 
by FLUKA and GEANT4. 


production data from HARP and E910 in Eq. |T|were sam¬ 
pled within their correlation errors. The sampling of these 
correlated parameters was performed via the Cholesky de¬ 
composition of the covariance matrix reported in Tab. 5 
of [ 23 ] ■ For each sampling of the c, coefficients, neutrinos 
were weighted with a factor 




da 

dp^dd-n 

da 

dp-n dO^ 


(g) 

w 


( 2 ) 


depending on the momentum (p^) and angle (0 T ) of their 
parent pion. c° are the best-fit values to the HARP-E910 
data-set. The FNR for different Ci varied within their co- 
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Fig. 12. The FNR results for the six configurations are shown as function of the neutrino energy. The effect of data-driven 
hadro-production uncertainties on the FNR computed using the Sanford-Wang parametrization for the six configurations is 
depicted. Histogram lines show typical individual-samplings of the d parameters of Eq. [l] Plain bullets correspond to the 
average value in each bin while the error bars represent the r.m.s. of the samplings. For each configuration, bottom plots (hollow 
bullets) show the ratios of the r.m.s. over the central values, providing an estimate of the fractional systematic error. 
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variance error-matrices are shown in Fig. [12] for the six 
considered configurations. 

For each configuration, in the top plots (plain bullets) 
the average value is shown with its error bar representing 
the r.m.s. of the samplings. Bottom plots (hollow bullets) 
show the ratio of the r.m.s. over the central value providing 
an estimate of the fractional systematic error. Uncertain¬ 
ties are rather large (5-7%) when considering the full area 
of the near detector at 110 m; they decrease significantly 
when restricting to the central region (configurations 2, 4 
and 6). In particular, in configuration 4, that is a realistic 
one from practical considerations, the uncertainty ranges 
from 2% at low neutrino energy and decreases below 0.5% 
-1.5% at neutrino energies above 1 GeV. The uncertainty 
is generally below 0.5% for a near site at 460 m. 

3.5 Conclusions for Section [3] 

Full simulations of the Booster beam were made anew, 
based on FLUKA and GEANT4 Monte Carlo simulations. 
Indications of the systematic error on the far-to-near ratio 
were obtained, showing characteristic behaviours. More¬ 
over, using the constraints from HARP-E910 data the un¬ 
certainties on the FNR, associated to hadro-production, 
were carefully estimated. 

Six configurations of the detector locations and sizes 
were considered. For a far detector on ground-surface and 
a near detector at the same off-axis angle the systema¬ 
tic error stays at 1-2% when the far and near transverse 
surfaces are matched in acceptance (configuration 4). Pro¬ 
vided the high available statistics and the large lever-arm 
for oscillation studies, the layout with baselines of 110 m 
and 710 m is considered in the following as the best choice. 


4 Detector Design Studies 

The location of the Near and Far sites is a fundamental 
issue in a search for sterile neutrino at SBL. Moreover the 
two detector systems at the two sites have to be as similar 
as possible. The NESSiE far and near spectrometer sys¬ 
tem were designed to match with a timely schedule and 
also exploit the experience acquired in the construction, 
assembling and maintenance of the OPERA spectrome¬ 
ters [28] that own an active transverse area of 8.75 x 8.00 

o 

m . 

The OPERA two large dipole iron magnets will be 
dismantled in 2015-2016 and possibly be re-used for the 
Vp disappearance study discussed in this paper. They are 
made of two vertical arms connected by a top and a bot¬ 
tom return yoke. Each arm is composed of 12 planes of 

5 cm thick iron slabs, interleaved by 11 planes of Resis¬ 
tive Plate Chambers (RPC) that provide the inner tracker. 
The magnetic field has opposite directions in the two arms 
and is uniform in the tracking region, with an intensity of 
1.53 T. 

Muons stopping in the spectrometers can be identified 
by their range. Their fraction can be maximized by in¬ 
creasing the depth of the magnets. This can be achieved 


by longitudinally coupling the two OPERA spectrometers, 
both at the far and the near sites, minimizing therefore the 
detector re-design. Their modularity allows to cut every 
single piece at 4/7 of its height, using the bottom part for 
the far site and the top part for near site. In this way any 
inaccuracy in geometry (the single 5 cm thick iron slab 
owns a precision of few mm) or any variation of the mate¬ 
rial properties with respect to the nominal ones (they are 
at the level of few percent) will be the same in the two 
detection sites. The near NESSiE spectrometer will thus 
be a clone of the far one, with identical thickness along 
the beam but a reduced transverse size. 

With the proposed setup a very large fraction of muons 
from CC neutrino interactions is stopped in the spectrom¬ 
eter. For this class of events the momentum is obtained 
by muon range. For higher energies the muon momentum 
is determined from the muon track sagitta measured in 
the bending plane. The charge of the muon can be deter¬ 
mined when its track crosses few RPC planes (> 3). The 
distributions of hit RPC planes are shown in Fig. [T3] for 
charged and neutral current events. 



Fig. 13. (color online) Number of crossed RPC planes for 
charged and neutral current events (5 cm slab geometry). 


In the positive-mode running of the Booster beam the 
antineutrino contamination is rather low (see Fig. |14[ > . In 
this case the use of the charge identification is limited 
and it can contribute only to keep the related systema¬ 
tic error under control and well below 1% since the mis- 
identification of the charge (mis—ID) of the Spectrome¬ 
ters is about 2.5% in the relevant momentum range (see 
the bottom plot of Fig. 24). The situation is quite dif¬ 
ferent for the negative-mode running where the neutrino 
contamination is rather high (see Section [6|. 

In the following the performances (efficiency and pu¬ 
rity) of the spectrometer with 5 cm thick iron-slabs are 
evaluated in terms of Neutral Current (NC) contamination 
and momentum resolution, and compared to a geometry 
with 2.5 cm thick iron slabs. 

With 5 cm slabs the fraction e of neutrino interactions 
in iron inducing at least one RPC hit (e = (> 1 RPC)/all) 
is 68%. The efficiency for CC and NC events is ecc = 86% 









A. Anokhina et al.: Search for sterile neutrinos in the disappearance mode at FNAL 


11 



E v (GeV) 


Fig. 14. (color online) The fluxes of the Booster neutrino beam 
(from [23]). 
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Fig. 15. CC efficiency (sec) and purity as a function of the 
minimum number of hit RPC planes for the two spectrometer 
geometries, 5 cm thick iron slabs (squares) and 2.5 cm thick 
slabs (circles). 


and £j vc = 20%, respectively. That corresponds to a frac¬ 
tion of NC interactions over the total number of interac¬ 
tion NC/all = 8.1%. With a minimal cut of 2 crossed 
RPC planes, the NC contamination is reduced to 4.2%; 
by requiring 3 RPC planes the NC contamination drops 
to 3.0%. 

Using 2.5 cm thick slabs, the fraction of neutrino events 
with at least one RPC hit increases for both NC and CC 
events. The CC efficiency and the NC contamination are 
both larger with respect to the 5 cm geometry. In Fig. [15] 
Ecc and the purity, CC/all = 1 — NC/all , are shown as a 
function of the minimum number of crossed RPC planes, 
for either slab thickness. 

At the same level of purity the efficiencies in the two 
geometries are comparable. No advantage in statistics is 
obtained with thinner slabs if the same NC contamination 
suppression is required. In conclusion the already available 
5 cm thick iron slabs can be adopted. It is worthwhile to 
note that purities and efficiencies have been extensively 
checked not to be spoiled by second-order systematic ef¬ 
fects due to the convolution of the fluxes and the neu¬ 
trino cross-section with the detector acceptances at the 
two sites. That keeps the systematic error due to the neu¬ 
trino detection below 1% while the most relevant contribu¬ 
tion to systematics remains the uncertainty on the fluxes]^] 

4.1 Track and momentum reconstruction 

The RPC digital read-out is provided on both vertical 
(Y) and horizontal ( X) coordinates using 2.6 cm pitch 
strips. Track reconstruction is made first in the two RPC 
projections (the XZ bending plane, and the YZ plane). 
Then, the two 2 D tracks are merged to reconstruct a 3-D 
event. 

For muons stopping inside the spectrometer the mo¬ 
mentum is obtained from the track range using the con- 
tinuous-slowing-down approximation 129] . The range dis¬ 
tribution is shown in Fig. |16| Similar conclusions can be 

3 More discussion is provided in the proposal [22] . 


drawn as for the number of crossed RPC planes (in Fig. 131, 
namely the CC efficiency and the reduction of the NC 
related background. For any muon track a parabolic fit 



Fig. 16. (color online) Muon range of reconstructed events, 
for CC plus NC (blue), CC (red) and NC (yellow) events. 


is performed in the bending plane to evaluate the track 
sagitta thus determining particle charge and momentum. 

In Fig. 17 (CC events) the reconstructed variables, 
namely the number of fired RPC planes and the range 
in iron, are plotted versus the muon momentum. A corre¬ 
lation is visible for both variables. The very strong muon 
momentum-range linear correlation allows to reach a sen¬ 
sitivity of few percent in the momentum estimation. 


4.2 Conclusions for Section 0] 

In the previous Section [3] it was shown that by adopting a 
realistic layout configuration, and by a proper choice of the 
fiducial volume at the near site, a reduction of the uncer¬ 
tainty on the Far/Near ratio to less than 2% level is pos- 
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Fig. 17. (color online) Number of crossed RPC planes (top) 
and range (bottom) versus muon momentum (CC events). 


sible, by plugging the data-driven knowledge on hadro- 
production (HARP and E910 [23] ). This is by far the dom¬ 
inant component. Possible other effects due to the running 
conditions of the detectors once installed can be kept un¬ 
der control (< 1%). It must be noted that the efficiency 
and acceptance of the spectrometers can be checked rou¬ 
tinely with good accuracy using the large amount of avail¬ 
able cosmic rays. Furthermore the detector in itself is very 
well mastered and understood due to its simplicity and to 
the extensive experience in running it underground on the 
CNGS beam. Furthermore, each original (from OPERA) 
iron slab will be used partly at the near and partly at the 
far site, providing the same geometrical and material com¬ 
position. That choice would provide a very constrained 
system at the two sites, with not only identical targets, 
but also similar geometrical frames and acceptances. The 
relative large statistical sample that could be obtained 
within configuration 4 would allow a careful control of the 
related systematic effects, by operating at different energy 
ranges, too. Finally, all the effects due to detector accep¬ 
tance and event-reconstruction have been evaluated to be 
within 1%. 


5 Physics Analyses and Performances 

The disappearance probability of muon-neutrinos, P{v^ —► 
nottyj, in presence of an additional sterile-state can be 
expressed in terms of the extended PMNS m mixing 
matrix (U a i with a = e,fi,T,s , and i = 1, ...,4). In 
this model, called “3+1”, the neutrino mass eigenstates 
vi, ... , 1/4 are labeled such that the first three states are 
mostly made of active flavour states and contribute to 
the “standard” three flavour oscillations with the squared 
mass differences Am 2 i ~ 7.5 x 10 -5 eV 2 and ~ 

2.4 x 10~ 3 eV 2 , where Am 2 ,- = mf — m 2 . The fourth mass 
eigenstate, which is mostly sterile, is assumed to be much 
heavier than the others, 0.1 eV 2 < Amh < 10 eV 2 . The 
opposite case in hierarchy, i.e. negative values of A 
produces a similar phenomenology from the oscillation 
point of view but is disfavored by cosmological results on 
the sum of neutrino masses [31] . 

In a Short-Baseline experiment the oscillation effects 
due to Z\m|i and Am 2 -, can be neglected since L/E ~ 
1 km/GeV. Therefore the oscillation probability depends 
only on Am 2 X and I / a 4 with a = e,n,T. In particular 
the survival probability of muon neutrinos is given by the 
effective two-flavour oscillation formula: 

^)l + bl = 1 - [4|tV| 2 (l - |tU 2 )] sin 2 A ™^ L , 

(3) 

where 4 11 /^ 41 2 (1 — 1 1 /^ 4 1 2 ) is the amplitude and, since the 
baseline L is fixed by the experiment location, the oscilla¬ 
tion phase is driven by the neutrino energy E. 

In contrast, appearance channels (i.e. v tl -+ v e ) are 
driven by terms that mix up the couplings between the 
initial and final flavour-states and the sterile state, yield¬ 
ing a more complex picture: 

Ve)% + B \ = 41 G/^ 4 1 2 1 C/ e 4 1 2 sin 2 (4) 

Similar formulas hold also assuming more sterile neutrinos 
(3 + n models). 

Since \U a i\ is expected to be small, the appearance 
channel is suppressed by two more powers in \U a i\ with 
respect to the disappearance one. Furthermore, since v e 
or appearance requires |C/ e4 | > 0 and > 0 , it 

should be naturally accompanied by non-zero v e and 
disappearances. In this sense the disappearance searches 
are essential for providing severe constraints on the theo¬ 
retical models (a more extensive discussion on this issue 
can be found e.g. in Section 2 of [32]). 

It should also be noticed that a good control of the 
v e contamination is important when using the v ^ —» v e 
for sterile neutrino searches at SBL. In fact the observed 
number of v e neutrinos would depend on the -+ v e ap¬ 
pearance and also on the v e —> v s disappearance. On the 
other hand, the amount of neutrinos would be affected 
by the -+ v s and v e —> transitions. However the 
latter term {v {1 appearance) would be much smaller than 
in the v e case since the v e contamination in v„ beams is 
usually at the percent level. In conclusion in the v jt disap¬ 
pearance channel the oscillation probabilities in either the 
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near or far detector are not affected by any interplay of dif¬ 
ferent flavours. Since both near and far detectors measure 
the same single disappearance transition, the probability 
amplitude is the same at both sites. 

Another important aspect of the analysis is related to 
the procedure of either rejecting or evaluating the pres¬ 
ence of a sterile component. The basic hypothesis of no- 
sterile oscillation (Ho) has been assumed against the pres¬ 
ence of something else (Hi). For Ho the analysis is sim¬ 
plified since the systematic errors on the FNR estimator 
are definitively under control when no-sterile component is 
included, as illustrated in the previous sections. In partic¬ 
ular the cross-section uncertainties, the hadro-production 
modeling, the beam-flux variations and their convolutions 
with the detectors’ acceptance have been checked: the sys¬ 
tematic error on FNR is 1 — 2%. Evaluation of p-values for 
Ho allows to set the possible presence of a sterile compo¬ 
nent. However, in order to estimate the power of an exper¬ 
iment exclusion plots should be also evaluated. The dis¬ 
tortion of FNR due to a sterile neutrino component with 
respect to the null hypothesis has been looked through, 
taking care of the correlations due to the systematic er¬ 
rors. In the following both procedures are depicted. 

The experiment sensitivity to the u ll: disappearance 
was evaluated by considering several estimators, related 
either to i) the muon produced in v < ^ c or to ii) the re¬ 
constructed neutrino energy. The muon momentum can 
be very effective when Hq hypothesis is checked to estab¬ 
lish the probability of the non-sterile component, i.e. the 
observation of something else. In such a case the simu¬ 
lation is limited to the standard processes and the FNR 
approach in the NESSiE environment is fully efficient. In¬ 
stead, to evaluate exclusion plots one needs to extract the 
oscillation parameters via Monte Carlo by looking at the 
FNR distortion in specific regions of the phase space. A 
new procedure based on the reconstructed muon momen¬ 
tum was also implemented to exclude regions defined by 
new “effective” variables. Using reconstructed measured 
quantities allows to keep systematic errors under control. 

In the second case (ii) the neutrino energy was recon¬ 
structed from 

E, u - ml/(2M) 

" 1 — (E^ — p^cos9)/M' 

valid in the Charge Current Quasi Elastic (CCQE) ap¬ 
proximation, M being the nucleon mass, E^ and the 
muon energy and momentum, respectively. 

We developed complex analyses to determine the sen¬ 
sitivity region that can be explored with an exposure of 
6.6 x 10 20 p.o.t., corresponding to 3 years of data col¬ 
lection on the FNAL-Booster beam. Our guidelines were 
the maximal extension at small values of the mixing angle 
parameter and the control of the systematic effects. 

The sensitivity of the experiment was evaluated per¬ 
forming three analyses that implement different techniques 
and approximations: 

• method I: a Feldman&Cousins technique (see Section 
V of [33]) with ad hoc systematic errors added to the 
muon momentum distribution; 


Table 4. Fiducial mass and baselines in configuration 4 for 
near and far detectors, used for the sensitivity analyses. 

Fiducial Mass (ton) Baseline (m) 

Near 297 FLO 

Far 693 710 


• method II: a Pearson’s y 2 test [29] with a full correla¬ 
tion matrix based on full Monte Carlo simulation and 
reconstruction; 

• method III: a new approach based on the profile likeli¬ 
hoods, often referred to as modified frequentist method 
or CL S 34], similar to that used in the Higgs boson 
discovery [35] . 

Throughout the analyses the detector configuration 
defined in Table 0] was considered. 

The distributions of events, either in E v or p^, normal¬ 
ized to the expected luminosity in 3 years of data taking 
(6.6 x 10 20 p.o.t.) with the FNAL-Booster beam running 
in positive focusing mode, are reported in Fig. [18| 

The study of the disappearance is reported in Sec¬ 
tion [6] with results obtained from method III. 


5.1 Sensitivity Analyses 

In the three analyses the two-flavour neutrino mixing in 
the approximation of one mass dominance was considered. 
The oscillation probability is given by: 


P&n ->^) = 1- sin 2 (26Q e ,„) sin 2 


( l-27Am 2 eu ,L[km] \ 

V E[ GeV] ) ’ 


where Z\m 2 etu is the mass splitting between a new heavy- 
neutrino mass-state and the heaviest among the three SM 
neutrinos, and 9 new is the corresponding effective mixing 
angle. 

In the selected procedures (Feldman&Cousins approach, 
X 2 test with Near-Far correlation matrix and CL S profile 
likelihoods) the evaluation of the sensitivity region to ster¬ 
ile neutrinos was computed at 95% C.L. Some recent use 
of more stringent Confidence Limits (even to 10 cr’s [36] ) 
was judged un-necessary, provided the correct and con¬ 
servative estimation of the systematic errors. Besides, we 
note that the measurement of muon tracks is a quite old 
and proven technique with respect to the more difficult 
detection and measurement of electron-neutrino interac¬ 
tions in Liquid-Argon systems. 


5.1.1 Method I (Feldman&Cousins technique) 

In method I the far-to-near ratio was written as Ri = 
Fi/(kNi), where F) and JV* are the number of events in the 
i-th bin of the muon-momentum distribution in the far 
and near detectors, respectively, and k is a bin-independent 
constant factor used to normalize each other the near and 
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Fig. 18. (color online) The total number of expected CC 
interactions seen by the Far detector at 710 m, as a function 
of E v (top) and p y (bottom), for the expected luminosity in 3 
years of data taking (6.6 x 10 20 p.o.t.) with the FNAL-Booster 
beam in positive-mode running. The sub-sample correspond¬ 
ing to the CCQE component is also shown. 


far distributions. For each value of the oscillation param¬ 
eters, sin 2 2 9 new and Z\m 2 etu , the y 2 is computed as 
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where i?o,i is the far-to-near ratio in absence of oscilla¬ 
tion and crfl 0 . is the quadratic sum of the statistical er¬ 
ror and a fixed, bin-to-bin uncorrelated, systematic error. 
In the Feldman&Cousins approach a 2\y' 2 ut (sin 2 29 new , 
Am 2 new ) = x 2 (sin 2 2 9 neun Am 2 new ) - xhin cut is applied. 
For every (sin 2 2d new , Am^ ew ) oscillated spectra were gen¬ 
erated and fitted to obtain the Xmin■ The distribution 
of Ax 2 (sin 2 20 new , Am\ ew ) is cut at 95% to define the 
(sin 2 20„ eu ,, Z\m 2 eu; ) exclusion region. The critical value 
on Ax 2 u t can be determined by either sampling the Ax 2 
distribution as in Feldman&Cousins or by applying the 
standard y 2 ut fixed-value for the 95% C.L. and two degrees- 
of-freedom. It was verified that in both cases the obtained 
results are very similar in the whole (sin 2 2 9 new , Am^ ew ) 
space. 


Results are shown in Fig. [Bl] for a set of ten simulated 
null experiments. In the top plot a systematic error e sys = 
0 was used. In the bottom plot a bin-to-bin uncorrelated 
systematic-error e sys = 0.01 was assumed (see Section 
12.1 of [ 52 ] for more details). 




Fig. 19. Top: sensitivity curves at 95% C.L. with ten simulated 
toy experiments and no systematic uncertainties. In the fit pro¬ 
cedure p y was used as observable with a cut of 500 MeV/c. 
10 4 sampling points, uniformly distributed in log scale, were 
generated. The expected integrated luminosity for 3 years of 
data taking (6.6 x 10 20 p.o.t.) with the FNAL-Booster beam 
in positive-mode running was assumed. Bottom: as above but 
using a bin-to-bin uncorrelated systematic-error e sys = 0.01. 


5.1.2 Method II (y 2 test with Near-Far correlation matrix) 

In method II the sensitivity to the disappearance was 
evaluated using two different observables, the muon range 
and the number of crossed RPC planes. The correlations 
between the data collected in the far and near detectors 
are taken into account through the covariant matrix of the 
observables. The y 2 is given by 

N N 

x 2 = EE ( kN i - F i) (M- 1 ).. (kNj - Fj ), (8) 

i=1 j—0 
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where M is the covariance matrix m of the uncertainties 
(statistical and bin-to-bin systematic correlations [58]). 

The Vp disappearance can be observed either by a 
deficit of events ( normalization ) or, also, by a distortion 
of the observable spectrum (s/iap<jff), which are affected by 
systematic uncertainties expressed by the normalization 
error-matrix and the shape error-matrix, respectively. The 
shape error-matrix represents a migration of events across 
the bins. In this case the uncertainties are associated with 
changes not affecting the total number of events. Conse¬ 
quently, a depletion of events in some region of the spec¬ 
trum should be compensated by an enhancement in oth¬ 
ers. Details of the model used for the shape error-matrix 
can be found in [ 22 ] . 

The distributions of the muon range and of the number 
of crossed planes were computed using GLoBES [39] with 
the smearing matrices obtained by the full Monte Carlo 
simulation described in Section |4] 

By applying the frequentist method the \ 2 statistic 
distribution was looked at in order to compute the sensi¬ 
tivity to oscillation parameters. Different cuts on the range 
and on the number of crossed planes were studied. Fur¬ 
thermore sensitivity plots were computed by introducing 
bin-to-bin correlated systematic uncertainties by consid¬ 
ering either 1 % correlated error in the normalization or 
alternatively 1 % correlated error in the spectrum shape. 

As a representative result the sensitivity computed us¬ 
ing the range as observable and taking the 1 % correlated 
error in the shape, is plotted in Fig. [20] Instead the nor¬ 
malization correlated-error would slightly reduce the sen¬ 
sitivity region around Z\m 2 euJ = 1 eV 2 . Moreover, the 
sensitivity region obtained by using the sum of CC and 
NC events is almost the same as that obtained with CC 
events only (see Section 12.2 in [22]). That proves that the 
result is not affected by the NC background events. 

5.1.3 Method III (CL S profile likelihoods) 

In the profile CL S method we introduce a new test-statistics 
that depends on a signal-strength variable. By looking at 
Eq. [fi] the factor sin 2 2 9 new acts as an amplification quan¬ 
tity for a fixed Am^ ew - Therefore a signal-strength /i can 
be identified with sin 2 2 6 new to construct the estimator 
function: 

r _ 1 ~ P ' s ^ n (1-27 Am new Lpnr/Ey) 

1 - M • sin 2 (1.27 Aml ew L Near /E v )' 

In a simplified way, for each Z\m 2 ew , a sensitivity limit 
on fi can be obtained from the p-value of the distribu¬ 
tion of the estimator / in Eq. [9] in the assumption of 
background-only hypothesis. 

That procedure does not correspond to computing the 
exclusion region of a signal, even if it provides confidence 
for it. The exclusion plot should be obtained by fully de¬ 
veloping the CL S procedure as described in Section 12.3 

4 Note that the shape analysis looks at the same distributions 
of Method I. 



Fig. 20. (color online) 95% C.L. sensitivity obtained using 
the muon range for all the interaction processes (Quasi Elas¬ 
tic, Resonant, Deep Inelastic Scatering): CC (black), CC+NC 
(red) and for CCQE events only (blue). 1% bin-to-bin corre¬ 
lated error in the shape is considered. The expected integrated 
luminosity for 3 years of data taking (6.6 x 10 20 p.o.t.) with the 
FNAL-Booster beam in positive-mode running was assumed. 

of 122] . However, since we are here mainly interested in ex¬ 
ploiting the sensitivity of the experiment, the procedure 
provides already insights into that. Its result comes fully 
compatible with the previous two analyses, which follow 
the usual neutrino analyses found in the literature. 

Moreover, following the same attitude, an even more 
aggressive procedure can be applied. Since the deconvo¬ 
lution from to E„ introduces a reduction of the in¬ 
formation, we investigated whether the more direct and 
measurable parameter, p can be a valuable one. In such 
a case Eq. [9] becomes: 

r = 1 - h ' sin 2 (1.27 Am 2 new L Far /p M ) 

1 - At' sin 2 (1.27 Am% ew L Near /p,j,) 

The corresponding sensitivity plot is shown in Fig. |21| 
It provides an “effective” sensitivity limit in the “effec¬ 
tive” variables Am 2 and the reconstructed muon momen¬ 
tum, p^rec- By applying the Monte Carlo deconvolution 
from p Mjr ec to E„ we checked that the “effective” Am 2 is 
simply scaled-off towards lower values, not affecting the 
mixing angle limil[^] The merits of the p^^ec are multi¬ 
ple: it is not affected as E v by the propagation error due 
to the deconvolution process, the systematics errors due 
to the reconstruction of the events (efficiency, acceptance, 
background) are directly included since it corresponds to 
a measured quantity, the estimator is powerful in identi¬ 
fying a possible new signal/anomaly. 

5.2 Conclusions for Section [5] 

The sensitivity curves obtained with different analyses 
prove the possibility to explore a very large region in the 

5 The scaled-off feature is evident from the comparisons 
of the sensitivity curves in the two cases, either E v or „c 
(Figs. 51 and 52 of the original proposal [22j). 
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Fig. 21. The sensitivity at 95% of C.L. obtained by computing 
the modified raster-scan method, in a CL S framework, and by 
using the reconstructed muon momentum as estimator (equa¬ 
tion (101). The expected integrated luminosity for 3 years of 
data taking (6.6 x 10 2 ° p.o.t.) with the FNAL-Booster beam 
in positive-mode running was assumed. A conservative cut of 
Pii,rec > 500 MeV/c was applied. 


mass-scale and mixing-angle plane, larger than other cur¬ 
rent proposals. Using a configuration with two (massive) 
detectors, one at 110 m on-axis, and one at 710 m off-axis 
(configuration 4, see Tabl e |4| , the achievable sensitivity 
curves are drawn in Fig. 22 for several C.L., compared 
to existing limits mm and the predicted sensitivities of 
the SBN project PH- A systematic error of 1% has been 
assumed and a conservative cut p^^ec > 500 MeV/c 

was applied. A sensitivity to mixing angles below 10 -2 
in sin 2 2 9 new can be obtained in a large region of Am 2 
around 1 eV 2 scale. 

It is noted that by applying more elaborated recon¬ 
struction algorithms than those used in the present anal¬ 
ysis the 500 MeV/c cut in p M>rec could be lowered to 
200 — 300 MeV/c. The exclusion region could then be sig¬ 
nificantly extended to Am 2 < 1 eV 2 , even if systematic 
errors should generally be larger and therefore detailed 
studies would be required to really access that limit. 

We demonstrated that a sophisticated statistical tool 
(method III) can be applied to get hints of new neutrino 
states at lower mass-scale than that achievable with the 
usual ones (methods I and II), by making use of a different 
estimator (the reconstructed muon momentum), less de¬ 
pendent of the Monte Carlo simulation. Thus we conclude 
that, on top of the exclusion limits, a robust confidence is 
accomplished on the identification of a possible new signal. 


6 Sensitivity for the antineutrino 
disappearance 


In negative-focusing mode the Booster beam contains a 
large neutrino component. In terms of flux the contami¬ 
nation amounts to 15% at 1 GeV, 30% at 1.5 GeV and 
surpasses the antineutrino flux above 2 GeV (see Fig. 23). 
In this energy region the measurement of the charge on 


Fig. 22. The sensitivity curves obtained using the modified 
raster-scan method (equation (|9|), in a CL S framework, for 
90%, 95% and 99.7% C.L., for an expected integrated luminos¬ 
ity of 3 years of data taking (6.6 x 10 2<) p.o.t.) with the FNAL- 
Booster beam in positive-mode running. An uncorrelated 1% 
systematic error and a conservative cut of p^.rec > 500 MeV/c 
were used. The filled area corresponds to the MINOS [8] result 
and the green curve to the MiniBooNE/SciBooNE [7] limit (at 
90% C.L.). The two dashed lines corresponds to the sensitivity 
predicted by the new SBN proposal [IS]) at 90% and 99.7% 
C.L. 


event-by-event basis is an efficient tool. Although a com¬ 
prehensive study of the spectrometers’ ultimate perfor¬ 
mance goes beyond the scope of this paper, we estimated 
the improvement on the final sensitivity by using their 
charge ID capability. We applied Method III of Section |5.1.3[ 
under some additional assumptions on the neutrino-anti¬ 
neutrino components. The contamination of the neutrino 
events resulting from the charge mis-ID probability (rj) 
was included in an uncorrelated way to the near and far 
data samples, with/without the assumption that only an¬ 
tineutrinos oscillate. The assumption that the oscillation 
phenomenon acts differently for neutrinos and antineu¬ 
trinos puts limits on their possible correlated oscillation 
probabilities. Though the resulting sensitivity curves cor¬ 
respond to the worst case scenario, they provide a decou¬ 
pled insight into the possible measurement of the antineu¬ 
trino disappearance at 1 eV mass-scale and small mixing 
angle, for the first time. 

For the antineutrino disappearance search two main 
results exist, an old one by the CCFR experiment [50] 
and more recent ones from the MiniBooNE [6] and Mini¬ 
BooNE/SciBooNE [41]. The MiniBooNE results come with 
a 20%-25% contamination of the intrinsic neutrino flux 
and were obtained assuming that only antineutrinos oscil¬ 
late while neutrinos do not. 

A total integrated luminosity corresponding to 3 years 
of running at the Booster in negative-focusing mode was 
considered. The number of events that could be collected 
at the far detector is displayed in Fig. [24] The neutrino 
sub-component is highly enhanced because of its larger 
cross-section with respect to the antineutrino one. 
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Fig. 23. (color online) The flux components of the Booster 
antineutrino beam (from I23j). 


To evaluate our sensitivity to the signal-strength esti¬ 
mator (Section 5.1.31 a set of four samples has been con¬ 
sidered, under different assumptions: 


1 ) the pure anti-neutrino sample (perfect rejection of neu¬ 
trino contamination); 

2 ) the anti-neutrino sample with the neutrino contamina¬ 
tion as determined by the spectrometers charge mis- 
ID. The same oscillation law is assumed for neutrinos 
and anti-neutrinos; 

3) as above, but assuming no oscillation for the neutrino 
contamination; 

4) the anti-neutrinos sample with full neutrino contami¬ 
nation (no charge ID), assuming no oscillation for the 
neutrino contamination. 


Fig. [25] shows the corresponding sensitivity curves, that 
have less and less power going from assumption (1) to 
(4), as expected. The most sensitive curve (case 1 and 
red line) corresponds to the pure sample of antineutrinos. 
The black line for case (2) is obtained by including the 
muon charge mis-ID in the collected event sample and in¬ 
creasing correspondingly the statistical errors associated 
to the neutrino and anti-neutrino components. The curve 
for case 3 (blue line) is obtained by assuming that the 
neutrino component, identified by the muon charge, do 
not oscillate, then decreasing the total sample while con¬ 
tributing to the statistical error. The purple line for (4) 
indicates the sensitivity in case the neutrino component is 
not identified and assumed not oscillating. 

The key feature of the charge identification is appar¬ 
ent since the quite small contamination coming from the 
mis-ID produces small corrections. In contrast lacking of 
charge measurement oblige to assume an oscillation pat¬ 
tern for the contamination component and reduces drasti¬ 
cally the amount of the equivalent statistical sample and 
the sensitivity region (case 4). 
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Fig. 24. (color online) The total number of foreseen z7 M CC 
interactions seen by the Far detector at 710 m, as function 
of the E v (top) and the p M (bottom), for the expected lumi¬ 
nosity in 3 years of data taking (6.6 x 10 20 p.o.t.) with the 
FNAL-Booster beam in negative-mode running. Separate an¬ 
tineutrino and neutrino data are shown, with also the sub¬ 
samples corresponding to the CCQE component (dashed his¬ 
tograms). In the bottom plot the mis-ID of the muon charge 
as provided by the spectrometers, is overlaid (black curve and 
left scale) [28]. The grey zone in the same plot corresponds to 
the expected CC neutrino contamination once measured the 
muon track (integrated 3200 over 58600 CC events in the 0.3 
- 5 GeV interval). 


7 Conclusions 


Existing anomalies in the neutrino sector may hint to the 
existence of one or more additional sterile neutrino states. 
A detailed study of the physics case was performed to 
set up a Short -Baseline experiment at the FNAL Booster 
neutrino beam exploiting the study of the muon-neutrino 
charged-current interactions. An independent measure¬ 
ment on complementary to the already proposed ex¬ 
periments on i/ e , is mandatory to either prove or reject the 
existence of sterile neutrinos, even in case of null result 
for v e . Moreover, very massive detectors are mandatory 
to collect a large number of events and therefore improve 
the disentangling of systematic effects. 
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Fig. 25. (color online) The four sensitivity curves (see text for 
explanation) obtained using the modified raster-scan method 
(equation |9j)), in a CL S framework, for 95% C.L., for an ex¬ 
pected integrated luminosity of 3 years of data taking (6.6 x 
10 2 ° p.o.t.) with the FNAL-Booster beam running in negative- 
focusing mode, r/ stays for the charge mis-ID probability. An 
uncorrelated 1% systematic error and a conservative cut of 
Pii,rec > 500 MeV/c were used. Previous results are also 

shown. The dashed green curve corresponds to the 90% C.L 
MiniBooNE [6] limit while the plain green curve corresponds 
to its improved analysis and data collection with also the Sci- 
BooNE data bh • 


The best option in terms of physics reach and fund¬ 
ing constraints is provided by two spectrometers based on 
dipoles iron magnets, at the Near and Far sites, located 
at 110 (on-axis) and 710 m (on surface, off-axis) from 
the FNAL-Booster neutrino source, respectively, possibly 
placed behind the proposed LAr detectors. 

A full re-use of the OPERA spectrometers, when dis¬ 
mantled, would be feasible. Each site at FNAL can host 
a part of the two coupled OPERA magnets, based on 
well know technology, allowing to realize “clone” detec¬ 
tors at the Near and Far sites. The spectrometers would 
be equipped with RPC detectors, already available, which 
have demonstrated their robustness and effectiveness. 

With that configuration one would succeed in keeping 
the systematic error at the level of 1 — 2% for the mea¬ 
surements of the interactions, i.e. the measurement of 
the muon -momentum at the percent level and the iden¬ 
tification of its charge on event by-event basis, extended 
to well below 1 GeV. 

The achieved sensitivity on the mixing angle between 
the standard neutrinos and a new state is well below 0.01 
for the v tJ mode. The measurement of the muon charge on 
event-by-event basis has been demonstrated to be very 
efficient for the estimation of possible disappearance an¬ 
tineutrino phenomena, for the first time at the level of few 
percents for the mixing angle and a mass scale around 1 
eV. 
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